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LONG-TERM  GOALS 

One  of  the  primary  research  goals  at  RSMAS  is  real-time  forecasting  of  both  Eulerian  fields,  such  as 
temperature  and  velocity,  and  Lagrangian  trajectories.  The  five  primary  components  are  (i)  MICOM, 
the  Miami  Isopycnal  Coordinate  Ocean  Model,  (ii)  satellite-derived  sea  surface  temperature  and  height 
fields  and  data  from  Lagrangian  drifters,  (iii)  an  Extended  Kalman  Filter  (EKE)  with  a  second-order 
Gauss-Markov  Random  Field  (GMRF)  model  for  spatial  covariances,  (iv)  a  random  flight  turbulence 
model  for  Lagrangian  trajectory  prediction,  and  (v)  contour-based  parameter  estimation  and  assimilation 
techniques. 

OBJECTIVES 

Documenting,  understanding,  and  predicting  ocean  variability  through  the  use  of  new  data  analysis  and 
assimilation  techniques. 

APPROACHES 

Our  data  analysis  and  assimilation  approaches  are  based  on  motion-compensated  space-time 
interpolation  algorithms,  state  space  reduction  techniques,  hodography,  and  multi-scale  field 
decomposition. 
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TASKS  COMPLETED 


Satellite-derived  SST  fields,  starting  from  1986,  are  being  assimilated  into  the  MICOM  using  the 
parameter  matrix  objective  analysis  algorithm.  Preliminary  results  indicate  a  more  realistic  Gulf  Stream, 
improved  estimates  of  both  the  Meridional  Overturning  Circulation  and  the  range  of  subtropical  mode 
water. 

Twin  experiments  using  our  MRF-based  EKF  for  assimilating  Sea  Surface  Height  (SSH)  into  a  multi¬ 
layer  MICOM  double  gyre  simulations  are  being  performed.  Prediction  errors  and  filter  performance 
were  evaluated  for  a  suite  of  updating  strategies,  subspace  projections/covariance  parameterizations, 
initial  conditions,  data  sampling,  and  model  variance  to  measurement  variance  ratios.  Experiments  with 
realistic  coastlines  and  topography  in  the  DAMEE-NAB  are  being  analyzed. 

An  algorithm  for  applied  Lagrangian  prediction  has  been  tested  in  model  simulations  (Ozgokmen  et  al., 
1998)  and  Gulf  Stream  search  and  rescue  simulations  with  real  drifter  data  (Schneider,  1998)  with 
encouraging  results  (with  A.  Griffa,  D.  Olson,  T.  Ozgokmen,  L.  Piterbarg  (USC)). 

Online  regional  and  global  data  sets  are  available  through  an  ftp  site  at  playin.rsmas.miami.edu.  In 
particular,  satellite-derived  quasi-global  sea  surface  temperature  fields  at  18  km,  2  day  resolution  for 
1985  to  mid  1998  are  now  online.  The  Gulf  Stream  Northern  Edge  Position  (GSNEP)  data  set  has 
been  extended  to  include  paths  from  1982-1996. 

Model-data  comparisons  for  sea  surface  temperature  (Wilson-Diaz  et  al.,  1998)  and  sea  surface  velocity 
(Han  et  al.,  1998)  have  been  completed  for  the  Arabian  Sea  and  Indian  Ocean,  respectively.  Model-data 
comparisons  studies  are  underway  with  Mike  Clancy  (FNMOC)  and  the  Navy  real-time  model. 

Semi-annual  to  interannual  variability  in  the  position  of  the  Gulf  Stream  was  investigated  as  a  function 
of  wind-  and  buoyancy- forcing  indicators  (Mariano  et  al.,  1998). 

RESULTS 

The  prognostic  fields  from  a  highly  nonlinear  multi-layer  MICOM  data  assimilation  experiments  exhibit 
stable  behavior  during  the  two-years  of  simulations  with  decreasing  forecast  errors  (manuscript  in 
preparation). 

Based  on  our  analysis  with  the  GSNEP  data  set  and  14  forcing  indicators,  we  find  that,  on  the  average, 
the  Gulf  Stream  moves  south  as  winds  increase,  as  temperature  decrease  and  the  large-scale  sea  surface 
thermal  gradient  increases.  The  south-to-north  displacement  from  the  first  PC  time  series,  that  explains 
46%  of  GSNEP  variability,  is  coherent  at  the  95%  confidence  level  with  (i)  the  root-mean-square  (rms) 
wind  stress  over  the  Gulf  Stream  on  time  scales  between  9  months  and  six  years,  (ii)  the  rms  wind  stress 
south  of  the  Gulf  Stream  on  5-6  year  time  scales,  and  (iii)  the  North  Atlantic  Oscillation  on  12-18 
months  time  scales.  The  first  PC  is  also  marginally  coherent,  80%,  with  the  Florida  Strait  cable  transport 
on  annual  time  scales  and  the  rms  wind  stress  south  and  north  of  the  stream  on  time  scales  greater  than 
one  and  two  years,  respectively  (Mariano  et  al.,  1998). 


The  meridional  component  of  the  second  PC  time  series,  that  explains  15%  of  GSNEP  variability,  is 
coherent,  at  the  95%  confidence  level  with  the  large-scale  thermal  gradient  across  the  Gulf  Stream  on  9- 
18  month  time  scales  and  (ii)  with  the  rms  wind  stress  south,  over,  and  north  of  Gulf  Stream  for  six 
months  to  six  year  time  scales.  The  second  PC  is  also  marginally  coherent,  80%,  with  the  Florida  Strait 
cable  transport  on  annual  time  scales.  The  first  two  PC  time  series  were  not  coherent  with  the  Southern 
Oscillation  Index  or  the  mean  latitude  of  the  large-scale  zero  wind-stress  curl  in  the  North  Atlantic.  The 
results  of  the  coherence  calculations  and  a  PCA  indicate  that  wind-forcing  and  buoyancy-forcing  are 
equally  important  for  determining  the  mean  latitude  of  the  Gulf  Stream  east  of  Cape  Hatteras  (Mariano 
et  al„  1998). 

The  first  mode  explains  64%  of  sea  temperature  variability  in  the  Arabian  Sea  and  is  associated  with 
wind- forcing.  The  second  mode  explains  31%  of  the  variability  and  is  associated  with  radiative  heating 
(Wilson-Diaz,  et  al.,  1998). 

Please  see  the  progress  report  by  Griffa  and  Ozgokmen  for  details  of  the  Lagrangian  work. 

IMPACT/APPLICATION 

Our  work  continues  to  be  well  reviewed  by  the  oceanographic  community  and  constitutes  significant 
advances  in  data  analysis  and  assimilation.  We  believe  that  this  work  will  lead  to  improved  ocean 
forecasts  that  will  save  fives  and  resources. 

TRANSITIONS 

The  lead  PI  presented  four  lectures  at  the  oceanography  division  of  NRL,  Stennis  Space  Center,  MS  on 
data  analysis  and  assimilation.  I  interacted  with  quite  a  number  of  researchers  during  my  week  long  visit. 
In  particular,  the  use  of  my  parameter  matrix  objective  analysis  for  motion-compensated  space-time 
interpolation  of  global  satellite  data  sets  was  discussed  with  Dr.  Barron  and  the  MRF-based  EKF  was 
discussed  with  Drs.  Shulman  and  Smedstad. 

RELATIONSHIPS 

Mariano  and  Chin  are  affiliated  with  the  RSMAS  remote  sensing  group  and  work  closely  with  the 
MICOM  modeling  group.  Strong  collaboration  with  RSMAS  scientists  A.  Griffa,  D.  Olson,  T. 
Ozgokmen,  as  well  as,  M.  Swenson  (AOMF),  and  F.  Piterbarg  (USC)  on  applied  Fagrangian  prediction 
will  continue  to  be  one  of  our  primary  near-term  research  activity.  We  are  working  with  A.  Moore 
(CS)  on  error  growth  and  contour-based  assimilation  and  with  F.  Shay  on  the  analysis  of  in-situ  ocean 
observations.  Model-data  comparison  studies  for  sea  surface  velocity  and  Fagrangian  prediction  studies 
are  underway  with  M.  Clancy  (FNMOC). 
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